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Abstract. Observational and theoretical results indicate 
that low-redshift BL Lacertae objects are the most likely 
extragalactic sources to be detectable at TeV energies. 
In this paper we present the results of observations of 4 
BL Lacertae objects (PKS0521-365, EXO 0423.4-0840, 
PKS2005-489 and PKS2316-423) made between 1993 
and 1996 with the CANGAROO 3.8 m imaging Cerenkov 
telescope. During the period of these observations the 
gamma-ray energy threshold of the 3.8m telescope was 
~2 TeV. Searches for steady long-term emission have 
been made, and, inspired by the TeV flares detected 
from Mkn 421 and Mkn 501, a search on a night-by-night 
timescale has also been performed for each source. Com- 
prehensive Monte Carlo simulations are used to estimate 
upper limits for both steady and short timescale emission. 

Key words: Gamma rays: observations - BL Lacs: indi- 
vidual: PKS0521-365 ,EXO 0423.4-0840, PKS2005-489, 
PKS2316-423 
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1. Introduction 

Only two extragalactic sources have been confirmed to 
emit 

gamma-rays above 300 GeV: Mkn 421 (Punch et al. 1992, 
Macomb et al. 1995, Petry et al. 1996) and Mkn 501 
(Quinn et al. 1996, Aharonian et al. 1997). In addition, 
an as yet unconfirmed detection has been made of 
1ES 2344+514 (Catanese et al. 1997a). All three have 
redshifts less than 0.05 and are, based on their values of 
log (Fx /F r ), classified as X-ray-selected BL Lacs (XBLs). 
Observations of more than 30 candidates of other source 
types over a range of redshifts have not yielded any other 
extragalactic TeV emitters (e.g. Kerrick et al. 1995). 
From theoretical considerations Stecker et al. (1996) have 
proposed that low redshift XBLs may be the only extra- 
galactic gamma-ray sources observable at TeV energies. 

The CANGAROO data set, containing data collected 
over more than 5 years, includes observations of four low 
redshift active galactic nuclei (AGN) of the BL Lacertae 
class (PKS0521-365, EXO 0423.4-0840, PKS2005-489 
and PKS2316— 423). In this paper we present the results of 
searches for TeV emission from these objects. A compar- 
ison of on-source and off-source regions of sky is used to 
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search for an excess of gamma-ray-like events from these 
sources over a typically two-week observing period. 

The detection of a large flare from the BL Lac Mkn 421 
by the Whipple Collaboration (Gaidos et al. 1996) has en- 
couraged us to analyze the data with particular empha- 
sis on shorter timescale flare searches. The Whipple re- 
sult has shown that at TeV energies, BL Lacs are ca- 
pable of extremely energetic flares on timescales of less 
than 1 day. Preliminary analysis of the Mkn 421 gamma- 
ray spectrum indicates that photon energies extend be- 
yond 5 TeV (McEnery et al. 1997). Similar strong, short 
duration flares at TeV energies have also been reported 
from Mkn 501, showing gamma-ray emission extending 
to energies of at least 10 TeV (Aharonian et al. 1997, 
Quinn et al. 1997). The lack of high energy turnover in 
the observed spectrum implies that the interaction of 
the gamma-rays with the cosmic infra-red background is 
at the lower end of expectations (cf Stecker et al. 1992). 
The results from Mkn 421 and Mkn 501 indicate that 
the CANGAROO 3.8 m telescope, with a relatively high 
gamma-ray energy threshold, is capable of detecting such 
extragalactic sources of TeV gamma-rays. 

2. The CANGAROO imaging atmospheric 
Cerenkov telescope 

The CANGAROO 3.8 m imaging telescope is located near 
Woomera, South Australia (longitude 137°47'E, latitude 
31°06'S, 160m a.s.l). The reflector is a single 3.8m diam- 
eter parabolic dish with f/1.0 optics. The imaging cam- 
era consists of a square-packed array of 10mm x 10mm 
Hamamatsu R2248 photomultiplicr tubes. The camera 
originally contained 224 photomultiplicr tubes, and was 
increased to 256 pixels (16 x 16 square array) in May 
1995. The tube centers are separated by 0.18°, giving 
a total field of view (side-side) of ~ 3.0°. The photo- 
cathode of each tube subtends 0.12° x 0.12°, giving a 
photo-cathode coverage of about 40% of the field of view. 
For a more detailed description of the 3.8 m telescope see 
Hara et al. 1993. 

In the current configuration an event trigger is gen- 
erated when a sufficiently large number of tubes (3~5) 
exceed their discriminator threshold. Individual tube dis- 
criminator levels are believed to be around 4 photo- 
electrons. Under these triggering conditions the current 
gamma-ray energy threshold of the 3.8 m telescope is es- 
timated to be ~1.5TeV. Prior to mirror re-coating in 
November 1996 (which includes all data presented in this 
paper) the energy threshold was somewhat higher. Us- 
ing Monte Carlo simulations we estimate that this energy 
threshold (as defined by the peak of the differential energy 
spectrum) was 2.5 TeV. When calculating integral fluxes 
and flux limits we define our threshold as 2 TeV. This fig- 
ure is obtained by re-binning the lower energy gamma-rays 
to produce an integral spectrum defined by a single power 
law with a sharp cutoff at the threshold energy. 



Since starting observations in 1991 the CANGA- 
ROO 3.8 m telescope has been used to observe a num- 
ber of galactic and extragalactic candidate TeV gamma- 
ray sources (Roberts et al. 1997). From these observations 
we have evidence for gamma-ray emission from three 
galactic sources — PSR1 706-44 (Kifune et al. 1995), the 
Crab nebula (Tanimori et al. 1994) and the Vela SNR 
(Yoshikoshi et al. 1997). 

As an indication of the sensitivity of the 3.8 m tele- 
scope to extragalactic sources we can use Monte Carlo 
simulations to estimate the response of the telescope to 
fluxes observed from Mkn 421 and Mkn 501. Assuming 
that the gamma-ray emission from Mkn 421 and Mkn 501 
extends up to 10 TeV the integral fluxes above 2 TeV are 
as follows. 

For Mkn 421: 
F(>2TcV) - 2.2 x 10~ 12 photons cm^s- 1 
(adopting the average 1996 Whipple flux with an assumed 
integral spectral index of —1.56, McEnery et al. 1997.) 

For Mkn 501: 
F(>2TeV) - 6.5 x 10~ 12 photons cm^s- 1 
(adopting the March 1997 HEGRA flux assuming an in- 
tegral photon spectral index of —1.49, Aharonian et al. 
1997.) 

If any of the sources examined in this paper were ca- 
pable of providing a steady flux at the level of Mkn 421 
or Mkn 501 it would be detectable by the CANGAROO 
3.8 m telescope, albeit at low significance. 

Observations by the Whipple group of Mkn 421 have 
shown that it is capable of extremely energetic flares 
on timescales of hours to days. A flare similar to 
that seen on 7 May 1996 by the Whipple telescope 
(McEnery et al. 1997) with a flux of F(>250GeV) - 5.7 x 
10~ 10 photons cm^s- 1 (F(>2TeV) - 3.8 x lO" 11 photons 
cm _2 s _1 for a —1.56 spectral index and maximum pho- 
ton energy of 10 TeV) lasting two hours would be easily 
detectable by CANGAROO at a significance of around 7<r. 

3. Data sample 

When selecting AGN as observation targets for the 3.8 m 
telescope we have used a number of criteria including the 
proximity of the source and measurements from X-ray 
and gamma-ray satellite experiments. As mentioned ear- 
lier, Stecker et al. 1996 have suggested that nearby XBLs 
are the most promising sources of detectable TeV gamma- 
ray emission. Multi-frequency studies of blazar emission 
(Sambruna et al. 1996) suggest that XBLs may be more 
compact and have higher electron energies and stronger 
magnetic fields than their radio selected counterparts. The 
current status of ground based TeV observations adds 
support to this belief. We present here the results of ob- 
servations on two XBLs and two radio selected BL Lacs 
(RBLs). While RBLs are perhaps less promising as TeV 
sources, well placed upper limits in the TeV range could 
help to confirm fundamental differences between XBLs 
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and RBLs. The four sources reported in this paper are 
described in the following sub-sections. 

3.1. PKS '052 1-365 

PKS0521-365 (RAj 2 ooo = 05 h 23 m 2.0, Decj 20 oo 
=— 36°23'2", z = 0.055) is a radio selected BL-Lac. It was 
first detected as a strong radio source more than 30 years 
ago (Bolton et al. 1964). The spectral energy distribution 
of PKS052 1—365 shows a double peaked structure that is 
typical of blazars (Pian et al. 1996). The first peak in the 
energy distribution occurs in the far IR (10 13 — 10 14 Hz) 
and is assumed to be from synchrotron emission from the 
electrons in the jet. The second component, possibly from 
IC emission from the same electron population, peaks at 
around lOOMeV. 

PKS052 1-365 was viewed by the EGRET experi- 
ment on the CGRO during 1992 from May 14 to June 4 
(Lin et al. 1995). A point source, consistent with the posi- 
tion of PKS0521— 365 was detected at a statistical signifi- 
cance of 4.3 a, with an integral source flux above lOOMeV 
of (1.8 ±0.5) x 10~ 7 photons cm _2 s _1 . The photon spec- 
trum from the EGRET observation can be fitted with a 
single power law 

dN/dE = (1.85±1.14)xl0- 8 (£yiGcV)- 2 - 16±0 - 36 pho- 
tons cm _2 s _1 GcV _1 

The hardness of the EGRET photon spectrum and the 
proximity of the source make PKS 0521-365 a candidate 
source for detectable levels of TeV gamma-ray emission. 

The CANGAROO data set for this object consists of 
52 on/off pairs of observations with a total of 89 hours of 
on-source and 84 hours of off-source data. 

3.2. EXO O423.4- 0840 

EXO 0423.4-0840 (RAj 200 o = 04 h 25 m 50.7, Dec. 72 ooo 
=— 08°33'43", z = 0.039) was reported as a serendipitous 
discovery as part of the high galactic latitude survey in 
the 0.05- 2.0keV range by EXOSAT (Giommi et al. 1991). 
Associating the source with the galaxy, MCG -01-12-005 
(incorrectly given in Giommi et al. 1991 as MCG +01 
12-005) and noting the high X-ray luminosity (> 10 43 
ergs s _1 ) of the source, Giommi et al. (1991) proposed the 
source as a candidate BL Lac object. This would make the 
source the third closest such object known (after Mkn 421 
and Mkn 501), and the closest in the southern hemi- 
sphere. We note however that Kirhakos & Steiner 1990 
determine a redshift of 0.0392 for the source they des- 
ignate IRAS 04235-0840B and identify with the HEAO 
source 1H 0422-086, and which they classify as a type 
2 Seyfert. Clearly higher resolution X-ray studies are re- 
quired to firmly establish the identity of this source. 



We have observed EXO 0423.4-0840 during October 
1996 obtaining a total raw data set comprising 20 hours 
of on-source and 17 hours of off-source data. 

3.3. PKS2005-489 

PKS2005-489 (RA 72 ooo = 20 /l 09 m 25.4, Dec 72 ooo = 
-48°49'54", z=0.071) is an XBL. X-ray measurements 
of PKS2005-489 by EXOSAT show extremely large flux 
variations on timescales of hours (Giommi et al. 1990). 
Initially reported as a marginal EGRET detection 
(von Montigny et al. 1995), a more accurate background 
estimation decreased the significance below the level 
required for inclusion in the Second EGRET Catalog 
(Thompson et al. 1995). However the fluxes above 1 GeV 
are more significant than those above 100 MeV suggesting 
that the source, in a part of the sky that has received rel- 
atively poor EGRET exposure, may indeed be detectable 
at higher energies (Lin et al. 1997). 

We have observed PKS2005— 489 during August of 
1993 and during August/September 1994 obtaining 41 
hours of on-source and 38 hours of off-source data. 

3.4. PKS2316-423 

PKS2316-423 (RA, 720 oo = 23 /l 19 m 05.8, Dec, 720 oo = 
-42°06'48", z= 0.055) is a radio selected BL Lac ob- 
ject. Assuming a magnetic field strength of B < 10~ 3 G, 
the emission from radio through to X-ray wavelengths is 
consistent with synchrotron radiation from electrons with 
E > 10 13 eV (Crawford & Fabian 1994). PKS2316-423 is 
not detected by the EGRET telescope on the CGRO. 

The CANGAROO 3.8 m telescope observed 
PKS2316-423 during July 1996 for a total of 26 hours 
of on-source data and 25 hours of off-source data. 

4. Image analysis and Monte Carlo simulations 

Prior to image analysis a data integrity check is performed 
to test for the presence of cloud or electronics problems. 
For a subset of observations tracking calibration is tested 
by monitoring changes in single-fold rates as local stars 
rotate through the field of view. Using this method the 
pointing direction can be inferred to an accuracy of around 
0.05°. 

For each event trigger the tubes associated with the 
image are selected using the following criteria: (i) The 
TDC of the tube must indicate that the tube has exceeded 
its discrimination threshold within 50 ns around the trig- 
ger time of the event, and (ii) The ADC signal in the tube 
must be at least 1 standard deviation above the RMS of 
background noise for that tube. 

An image is considered suitable for parameterization 
if it contains more than 4 selected tubes, and if the total 
signal for all tubes in the image exceeds 200 ADC counts 
(around 20 p.e.). About 25% of raw images are rejected 
by these two selection conditions. Surviving images are 
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parameterized after Hillas 1985, with the gamma-ray do- 5. Results 
mains for our observations being : 



0.5° < Distance < 1.1° 
0.01° < Width < 0.08° 
0.1° < Length < 0.4° 
0.4 < Concentration < 0.9 
alpha < 10° 

The cumulative percentages of events passing each se- 
lection condition are given in table 1. The numbers shown 
are based on image cuts applied to Monte Carlo gamma- 
rays and protons. The efficiency of the cuts for the Monte 
Carlo protons is consistent with that seen for real data. 



Cut 


Gammas 


Protons 


Raw data 


100% 


100% 


ADC sum>200 


90.0 


77.5 


Distance 


65.4 


44.6 


Width 


54.8 


20.4 


Length 


54.0 


16.5 


Concentration 


40.3 


10.9 


alpha 


38.8 


1.3 



Table 1. Cumulative percentages of events passing gamma-ray 
selection cuts. 



To estimate fluxes and upper limits for our data we use 
an exposure calculation based on a detailed Monte Carlo 
simulation of the response of our telescope to gamma- 
ray initiated EAS. Simulations have been based on the 
MOCCA simulation package (Hillas 1982) which models 
all relevant particle production processes and includes at- 
mospheric absorption effects for the Cerenkov photons 
that arc produced. The energies of simulated gamma-ray 
primaries were selected from a power law with integral ex- 
ponent — 1.4, with a minimum primary energy of 500 GeV. 
Core distances were selected from an appropriate distri- 
bution out to a limiting core distance of 250 m. Simulated 
gamma-ray images were then subjected to the same se- 
lection criteria as the real data. Fluxes and upper limits 
are calculated by comparing measured gamma-ray excess 
rates to those predicted by the simulations. If we assume 
that our telescope model is correct, this method of flux 
estimation has only two free parameters — the source 
spectral index and source cutoff energy. For total flux cal- 
culations the estimate of cutoff energy is not critical - 
the generally steep nature of gamma-ray source spectral 
indices ensures that the bulk of the flux is around the 
threshold energy. 
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Fig. 1. Scatter plot of night by night 2a flux limits for each of 
the sources (indicated by crosses). Also shown are the 2a flux 
limits for the total data set for each source (stars). The two 
dashed lines show the >2TeV fluxes from Mkn 421 (left) and 
Mkn 501 (right) (see text for details). 



The total data set for each source has been tested 
for the presence of gamma-ray signals. The significance of 
gamma-ray excesses have been calculated using a method 
based on that of Li & Ma 1983: 



+ N f fin 



Non+Noff 

Noff 



N on + Noff 



1/2 



(1) 



where S is the statistical significance and (3 is the ratio 
of events in the on-source observation to those in the off- 
source observation in the range 30° < alpha < 90° (where 
alpha is the image parameter describing image orienta- 
tion) . The values of N on and N Q f / are calculated from the 
gamma-ray selected (alpha < 10°) data for the on and 
off-source observations respectively. This method tends to 
slightly overestimate the significances because it does not 
account for the statistical uncertainty in calculating the 
value of (3. Based on Monte Carlo simulations we esti- 
mate that this effect is small — and less than the typical 
systematic differences caused by the variations in param- 
eter distributions between on and off-source observations. 
In the total data set for each source no significant excess 
is seen for those events in the gamma-ray domain. The 
calculated excesses are -0.27er (PKS0521-365), -0.99a 
(EXO 0423.4-0840), -0.91er (PKS2005-489) and +0.22cr 
(PKS2316— 423). Upper limits to steady emission have 
been calculated after Protheroe 1984 and are shown in 
the scatter plot in Fig. 1. 



Roberts et al. 



5 



PKS0521 


PKS2005 


PKS2316 


£!jyvvJU4Zo 


MJD 


F 


MJD 


p 


MJD 


F 


ivi j u r 


yooo. r 


9 Q 

z.y 


noon c; 

yzzu.o 


A 

4.0 


1U0U4.0 


4.o 


lUotw./ z.l 


Qftftl 7 

yooo. t 


ft ^ 

D.O 


Q991 ft 


O.O 


lUoUO. / 


9 7 
Z. ( 


lUoOo. / l.o 


Qft&A ft 

yoo4.o 


ft 7 


0.999 
yzzz.o 


7 1 


1 Ci^C\ft 7 
lUoUO. / 


O.l 


lUoby. / o.l 


yooo. < 


Q Q 

o.y 


yo i D. / 


A Q 

4.y 


lUou ( .o 


9 /I 
Z.4 


Woil). ( 0.4 


QfiQfi 7 

yooo. ( 


O.l 


QC77 7 

yo ( ( • i 


A 8 
4.o 


lUoUo. / 


1 £ 
1 .o 




yoo / . < 


A Q 
4.o 


ncn7 c 

yOy < .0 


7 n 


lUoll.O 


O.O 




yooo. ( 


Q 7 
O. / 


ncnn p; 

yoyy.o 


9 

Z.D 








yoyu. t 


^.O 


youu.o 


O.O 








yoy±. t 


7 1 


youi.o 


q 
o.y 








Q71 *3 

y / ±o.o 


A 0, 


youz.o 


ft 7 
O. / 








9714.7 


10.0 


9604.6 


8.9 








9715.6 


7.2 


9605.6 


4.1 








9717.6 


4.8 












9718.6 


8.2 












9719.6 


3.0 












9720.5 


4.2 












9722.6 


6.6 













Table 2. The 2cr flux upper limits, F(>2TeV) (xlCT 12 photons 
cm -2 s" 1 ) and approximate MJD (-40000) for each on-source 
observation. 



We have also searched our data set for gamma-ray 
emission on a night by night basis. In general our observa- 
tions of a source consist of a long (several hours) on-source 
run, with a similar length off-source run, offset in RA to 
provide the same coverage of azimuth and zenith. The flare 
search has been performed by calculating the on-source 
excess for each pair of on/off observations each night. In 
cases where there is no matching off-source run, an equiv- 
alent off-source run from another nearby night is used. 
Figure 2 shows the distribution of on-source significances 
for all three sources. There is no evidence for gamma-ray 
flares on the timescale of <~ 1 night for any of the sources. 
The most significant nightly excess (from PKS0521-365) 
has a nominal significance of 3.7<r but after allowing for the 
number of searches performed this significance is reduced 
to less than 3cr. The upper limits to gamma-ray emission 
for these observations are shown in Fig. 1. We have also 
included, in table 2, a list of upper limits to gamma-ray 
emission for each individual observation. 



6. Discussion 

The interpretation of upper limits from BL Lacs is dif- 
ficult, because at the present time there is no detailed 
model to predict TeV fluxes from the different classes 
of BL Lacs. Stecker et al. 1996 have attempted to predict 
TeV fluxes from a number of nearby XBLs by assuming 
that all XBLs have similar spectral characteristics to the 
known TeV emitter Mkn 421. Macomb et al. 1995 argue 
that the observed X-ray/TeV gamma-ray flux increases 




Fig. 2. Distribution of the significances of night by night ex- 
cesses for all sources. 



from Mkn 421 during flares indicate that the TeV gamma- 
rays are produced by a synchrotron self-Compton (SSC) 
mechanism. Stecker et al. 1996 note that the Compton 
emission from the electrons in the jet has a similar spec- 
trum to the synchrotron component, but upshifted by 
the square of the maximum Lorentz factor (estimated to 
be ~ 10 4 - 5 for Mkn 421) of the electrons. For Mkn 421 
and PKS2155-304 (both XBLs detected by EGRET) 
the luminosities in the Compton and synchrotron regions 
are nearly equal. Assuming that other XBLs also have 
Lc/L S y n <~ 1, Stecker et al. 1996 derive the following re- 
lationship between X-ray and TeV source fluxes : 

VTevFTeV ~ v x F x 

where v is the frequency and F the flux in each en- 
ergy band. Using an assumed E~ 2 - 2 spectral index for all 
XBLs, and absorption of TeV gamma-rays in the inter- 
galactic infra-red based on an average of Models 1 and 
2 from Stecker & de Jager 1997, they estimate TeV fluxes 
above 1 TeV from a number of nearby XBLs based on EX- 
OSAT X-ray flux measurements. It should be noted that 
a more recent paper (Stecker & de Jager 1998) concludes 
that the absorption of TeV gamma-rays in the infra-red 
background is overestimated in Stecker & de Jager 1997. 
Allowing for this, the predicted flux for PKS2005-489 
above 1 TeV is F(> ITeV) = 1.3xl(T 12 photons cm^s" 1 . 
Above 2 TeV (the threshold of CANGAROO in this anal- 
ysis) this flux would be F(> 2TeV) = 0.34 x 1(T 12 pho- 
tons cm _2 s _1 , below the flux sensitivity of the measure- 
ment made in this paper. The photon spectral index as- 
sumed in Stecker et al. 1996 is now incompatible with re- 
cent measurements of Mkn 421 and Mkn 501. It is also 
possible that the SSC mechanism is not primarily respon- 
sible for the TeV gamma-ray emission and a number of 
other mechanisms have been suggested (Mannheim 1993, 
Dermer & Schlickeiser 1994 and references therein). 
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Of the other predictions of Stecker et al. 1996, 
Kcrrick et al. 1995 derive an upper limit to the >0.3 TeV 
emission for the XBL 1ES 1727+502 that is above the cal- 
culated flux. It is also worth noting that a deep exposure 
of the RBL BL Lacertae with the Whipple telescope did 
not yield a detection (Catanese et al. 1997b). 

Current and future observations of a range of BL Lacs 
by ground based Cerenkov telescopes should help to clarify 
the mechanisms for the production of high energy photons 
in these sources. 

7. Conclusion 

Analysis of CANGAROO data shows no evidence for long- 
term or short-term emission of gamma-rays above 2 TeV 
from the BL Lacs PKS0521-365, EXO 0423.4-0840, 
PKS2005-489, PKS2316-423. The 2a upper limits to 
steady emission are 1.0 x lO^cm"^" 1 (PKS0521-365), 
1.1 x 10" 12 cm- 2 s- 1 (EXO 0423.4-0840), 1.1 x 
10" 12 cm- 2 s- 1 (PKS2005-489), and 1.2 x lO^cm"^" 1 
(PKS2316-423). For the XBL PKS2005-489 the flux 
limits presented in this paper do not constrain the 
TeV emission levels predicted by the simple model of 
Stecker et al. 1996. 
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